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Abstract 

The Finite Fourier Transformation matrix (F.F.T.) plays a central 
role in the formulation of quantum mechanics in a finite dimensional 
space studied by the author over the past couple of decades. An out- 
standing problem which still remains open is to find a complete basis 
for F.F.T. In this paper we suggest a simple algorithm to find the 
eigenvectors of F.F.T. 


Talk presented in the Second International Workshop on Harmonic Oscillators 
held in Cocoyoc , Mexico during 23-25 March, 1994, to appear in the 
Proceedings . 


■'i * * 


239 



I. INTRODUCTION 


The finite Fourier transform matrix (F.F.T.) plays a fundamental 
role in many contexts and has been studied extensively [1-3]. It is 
central in the discussions on finite dimensional quantum mechanics based 
on Weyl's commutation relations 14] studied by the author in a series of 
publications [5] . The eigenvalues of this matrix were determined by 
Schur [1] and a simple argument to recover this result has been given 
earlier [6], The calculation of the eigenvectors is not straight- 
forward and many methods have been given in particular, by Mehta [7]. 

In Section IV, we present a new algorithm to find the eigenvectors. 

II. EIGENVALUES OF S 

The F.F.T. matrix S, which is unitary, is defined by 

_ 1 r 2TFi 

S ~ ~ exp L a 6 ] , 

ag /n r n 

a, 3 = 0,1,2, . . . n-1 (2.1) 


i * /-I 


and has many interesting properties 

1) (s 2 ) . = I» - 6 q 

ag aB a + g, o 

(mod n) 


( 2 . 2 ) 


2 2 
Since S f = f , for a vector f with n components, S is 

a -a mod n, a r > 

called the parity operator 


2) (S ) Q = 6 Q 

aB aB 


(2.3) 

like the usual Fourier transform. 

3) The matrix S, which is by definition a symmetric matrix will 
diagonalize any circulant matrix. 
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From Equation (?-3), it is clear that the eigenvalues of S are simply 


±1 and ±i. There is then a degeneracy of the eigenvalues. The first 
problem will be to determine this. Luckily, Equations (2.1)-(2.3) can be 
repeatedly used to fix this [6] . If k^ s k^ > k^and k^ denote the multi- 
plicity of the eigenvalues taken in the order (1, -1, i, -i), Equation 

(2.1) implies that 

- n-1 ^ . 0 2 

1 ^ r 1 

Tr S = — I [ exp ] 

rv l=o n 


= | (1 + i) [ 1 + exp ( -yli ) ], (2.4) 

and hence 

Tr S = - k 2 > + i(k 3 - k^) 

= 1 for n = 4k + 1 , 

= 0 for n = 4k + 2 , 

= i for n * 4k + 3, 

= (1 + i) for n = 4k, 

k = 0,1,2, . . . 0.5) 


From Equation (2) we infer that 

Tr S 2 = (k l + k 2 ) - (k 3 + k 4 ) 

- 1 for n odd , 

= 2 for n even. 

We also have 

Tr = n = k^ + k 2 + k^ + k^. 

Equations (2.5), (2.6) and (2.7) can be used to solve for k^ k 2 , k 3 and 

k. and one finds that 
4 
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n - 4k + 1 


n = 4k + 2 


n = 4k + 3 


n = 4k 


k k + 1 


k, k 

4 

III- EIGENVECTORS OF S 


k + 1 
k + 1 
k 
k 


k + 1 
k + 1 
k + 1 
k 


k + 1 
k 
k 

k - 1 


Let us decompose S into its primitive 


S 


4 i 
E i J B(j), 

j-1 


idempotents as 


where 


B( 1 ) 

= }s + 

J <1 - I') 

B(2 ) 

= -}c + 

1 <1 + I') 

B(3) 

= -{s + 

i a - !•) 

B(4) 

= + 

4 (I + 


(3.2) 


1 2tt 

C = — cos ( — a8 ) 

ag n 

1 . , 2tt . 

s a6 ' H sl " ( ~ “ 6 >■ 

0.3 = 0,1,2, ...n-l (3.3) 

It is easily verified that 

S B(j) = i j B(j), (3.4) 

thus the nonzero columnus of B(j) yield the eigenvectors of S with eigen- 
value i^ . Also, ill analogy with the standard case, Mehta [7] hap been 
able to express these eigenvectors in terms of Hermite functions with 
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discrete arguments. 


IV. EIGENVECTORS OF S ; AN ALTERNATE METHOD 


Since the 

F.F.T. matrix S 

satisfies Equation (2.1) 

we construct 

the matrix [ 10 ] 






T = 

3 2 

s + s S J 

d 

+ 

s s 2 + 

d d 



= 

I' (S + 

S d ) 

+ (S + S d ) 


(4.1) 

where 






s d = 

diagonal 

s . 



(4.2) 

We find that 






S T 

= s (s 3 + 

s 2 s 

2 3 

. + S S . + S J 
d d d 




= (I + s 

3 s d 

2 2 

+ s s, + 

d 

s s’) 



= (S 3 d + 

s 3 

+ s 2 s_, + s 

d 

s d> s d 



= T S d . 




(4.3) 


If T is nonsingular, 

T + S T = S_, (4.4) 

d 


Therefore, the columns of T automatically provide the eigenvectors of 
S. The degenerate eigenvectors of S corresponding to the repeated eigen- 
values can be made orthonormal by using Gram-Schmidt process. This will 
render T unitary. While the process is quite general, we shall ; illustrate 
this for some special cases 
case of in = 2 


S 



(4.5) 
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and 



2 2 

Since S = S = I, 
a 

We get from Equation (4.1) 

T = 2 (S + S d ) , 



(4.6) 

(4.7) 


(4.8) 


We unitarized matrix of the eigenvectors of S is therefore 


+ 1 


U 2 = 


case of n = 3 


S = 


^ I'/i (SI + i) 


/3 


e = exp 


2tt i 


From Equation (2.{j) we see that 
/ 10 0 



1 - (/2 + 1 ) 




(4.9) 


(4.10) 


(4.11) 


one finds from Equation (4.1) that the unitarized matrix of the eigen- 
vectors of S is 
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/2 Si (Si + 1 ) 


Si + 1 


1 1 + {L i^3 + /3 

Si 

1 1 ± ^ -1/3 + /3 

Si 


(4.12) 


In this case we have 


s = I 


0 0 


0 0 


0 -1 


0 0 


It is easily calculated that 


-1 2i 


1 -1 -1 


-1 -2i 


The first two column vectors correspond to the eigenvalue +1, the third 


to -1 and the last to -i. 
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By a simple use of Gram-Schmidt orthogonalization procedure one 


can find the unitarized matrix corresponding to the eigenvectors of S 


as 


U, 


1 

/2 Jk (/4 + 1 ) 




3 

0 


0 

1 



i/6 

1 

-2/2 

-/3 

0 

1 



-ifi 


\ 

/ 


.16) 
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